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ABSTRACT 

We describe the infrared properties of sources detected over ~36 deg 2 of sky in the GAMA 15-hr 
equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H- 
ATLAS) and Wide-field Infrared Survey ( WISE). With 5cr point-source depths of 34 and 0.048 mJy 
at 250 /tm and 3.4 /im, respectively, we are able to identify 50.6% of the H-ATLAS sources in the 
WISE survey, corresponding to a surface density of ~ 630 deg~ 2 . Approximately two-thirds of these 
sources have measured spectroscopic or optical/ near-IR photometric redshifts of z < 1. Fbr sources 
with spectroscopic redshifts at z < 0.3, we find a linear correlation between the infrared luminosity 
at 3.4 /mi and that at 250 /Am, with ±50% scatter over ~ 1.5 orders of magnitude in luminosity, 

~10 9 - 10 10 * 5 L 0 . By contrast, the matched sources without previously measured redshifts (r >20.5) 
have 250-350 fim flux density ratios that suggest either high-redshift galaxies ( z > 1.5) or optically 
faint low-redshift galaxies with unusually low temperatures (T < 20). Their small 3.4-250 /Ain flux 
ratios favor a high-redshift galaxy population, as only the most actively star- forming galaxies at low 
redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large 
AGN fraction (~ 30%) in a 12 /Am flux-limited subsample of H-ATLAS sources, also consistent with 
there being a significant population of high-redshift sources in the no-redshift sample. 

Subject headings: surveys - cosmology: observations — galaxies: high-redshift - infrared radiation - 
galaxies: statistics - galaxies: general 
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1. INTRODUCTION 

The advent of submillimeter astronomy has opened a 
new window into the Universe, allowing us to probe dusty 
star- forming galaxies at high redshift, as well as cold 
dust in nearby galaxies. Along with the recent success 
of the balloon-borne BLAST experiment (Pascale et ah 
2008; Eales et al. 2009), the launch of the Herschel Space 
Observatory 1 (Pilbratt et ah 2010) allows us to probe 
wavelengths from 55 to 672 /xm to a spatial resolution 
of < 10". Its largest open-time key project, the Herschel 
Astrophysical Terahertz Large- Area Survey (H- ATLAS, 
Eales et al. 2010a), will observe ~550 deg 2 of sky and 
detect more than 300,000 galaxies. 

Based upon a preliminary cross-identification with 
the Galaxy and Mass Assembly survey (GAMA, 
Driver et al. 2011) and Sloan Digital Sky Survey (SDSS, 
Abazajian et al. 2009), ~5G% of the sources detected in 
the H- ATLAS survey are at z < 1 (Smith et al. 2011). In 
addition, an analysis of the FIR colors of the remaining 
H-ATLAS so luces with flux densities above 35 mjy at 
350 pu i and > 3cr detections at 250 fim and 500 pm sug- 
gests an average redshift of z ~ 2 (Amblard et al. 2010; 
Lapi et al. 2011). 

Another method of constraining the redshift distri- 
bution of FIR-selected sources is to look for counter- 
parts in the near and mid infrared (NIR and MIR). 
Fleuren et al. (in prep) have performed source match- 
ing to survey data from the VISTA Kilo-degree Infrared 
Galaxy survey (VIKING, Sutherland et al., in prep), but 
these data only extend to ~ 2 fim. Starting in December 
2009, the Wide-field Infrared Survey Explorer 2 (WISE, 
Wright et al. 2010) began its mission to observe the en- 
tire sky in four bands, ranging from 3.4 to 22 pm, at 
< 12" resolution. Following the first public data release 
in April 2011, ~ 24,000 deg 2 of IR images and source cat- 
alogs became available to the public, including 36 deg 2 
of sky in the equatorial plane covered by the H-ATLAS 
survey. 

The purpose of this letter is to describe the in- 
frared properties of identified WISE counterparts to H- 
ATLAS sources within the GAMA 15- hr field (G15). 
Throughout we will assume a concordance cosmology 
with Ho = 71 km s -1 Mpc -1 , Q, m = 0.27, and Ha = 0.73 
(Spergel et al. 2007). In a subsequent paper, we will re- 
port WISE cross-identifications over the entire H-ATLAS 
area and perform spectral energy distribution (SED) fits 
to the matched sources. 

2. DATA AND METHODOLOGY 

Herschel observations of the G15 field include imag- 
ing data at 250, 350, and 500 pm from the SPIRE in- 
strument, with respective beam FWHM of 18"1, 25'/2, 
and 36" 6 (Griffin et al. 2010), as well as 100 and 160 pm 
imaging from the PACS instrument (Poglitsch et al. 
2010). The field subtends approximately 12° in right 
ascension and' 3° in declination and is centered on the 
equatorial plane at a = 14 h 30 m. The H-ATLAS G15 
source catalogue (Dunne et al., in prep) is constructed 

1 Herschel is an ESA space observatory with science instruments 
provided by European- led Principal Investigator consortia and with 
important participation from NASA. 

2 For a description of the WISE mission, see 

http: // wise2.ipac.caltech.edu /docs / release /prelim/ expsup/ 


using the same method as the Science Demonstration 
Phase catalogue (Rigby et al. 2011). It contains 27,481 
sources detected at > 5cr in any of the three SPIRE bands 
and reaches point-source depths of 34, 40, and 44 mJy 
at 250, 350, and 500 pm, respectively. 

The WISE first public data release contains four-band 
coverage of the entire G15 field to 5 a point-source depths 
of 0.048, 0.10, 0.73, and 5.9 mJy at 3.4, 4.6, 12, and 
22 pm. The angular resolution in these bands is 6"1, 6"4, 
6"5, and 12?0, respectively (Wright et al. 2010). There 
are ~ 240,000 WISE sources within G15 detected at > 7a 
in at least one of the four WISE bands, corresponding to 
~0.14 WISE sources per 10"-radius aperture. In the 
WISE preliminary release catalog, close pairs of sources 
are not deblended for separations < 9", so when a back- 
ground or foreground source appears near an H-ATLAS 
source position, it will often be blended with the true 
counterpart. 

In addition to the infrared data from WISE and H- 
ATLAS, the G15 region has spectroscopic redshifts from 
the GAMA survey (r < 19.4, Driver et al. 2011), and 
photometric redshifts obtained using optical/NIR pho- 
tometry from SDSS, VIKING, and the UK Infrared 
Deep Sky Survey Large Area Survey (UKIDSS-LAS, 
Lawrence et al. 2007). These photometric redshifts are 
derived following Smith et al. (2011) and have typical 
redshift uncertainties of a z /z ~ 0.15. Of the 27,481 H- 
ATLAS sources in the region, 16% have reliable spectro- 
scopic redshifts and 39% have photometric redshifts. 

As a result of the high space density of 3.4 pm sources, 
a naive matching to the H-ATLAS source positions will 
result in a non-negligible number of misidentifications. 
A significant fraction of these misidentifications will be 
foreground stars that are blue in all WISE bands and 
undetectable in H-ATLAS, so we can reduce the con- 
tamination rate of our matched catalog by considering 
only WISE sources with [3.4] - [4.6] > 0. The magnitude 
distribution of all background/ foreground WISE sources 
is such that ~ 20% are bluer than this limit, as compared 
to < 1% of sources within 10" of an H-ATLAS position. 

To estimate the fraction of H-ATLAS sources with 
detectable counterparts in the remainder of the WISE 
source catalog, we use the method of Fleuren et al., in 
prep, where the detection rate is given by, 

Qo = i-f, (l) 

where S is the fraction _of unmatched positions in the 
H-ATLAS catalog and B is the fraction of unmatched 
random positions. Matching all WISE sources within 
10" (within which we expect > 99% of the true matches 
to lie) and using 10 5 random positions within the G15 
field, we find Q 0 = 0.632 =b 0.004. This is larger than 
the value found in the SDSS galaxy catalog (Qo = 0.583, 
Smith et al. 2011), but smaller than in the VIKING Ks 
band (Q 0 = 0.75, Fleuren et al., in prep). By contrast, 
we find Qo = 0.012 ± 0.002 within the subset of WISE 
sources with [3.4] - [4.6] < 0, suggesting that this color 
cut was effective in removing stars. 

In Fig. 1, we use Eq. 1 to estimate the IR detection rate 
as a function of FIR flux density. For sources brighter 
than 100 mJy at 250 pm, we find a WISE 3.4 pm coun- 
terpart > 95% of the time. The majority of these sources 
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Fig. 1. — NIR and MIR detection rates plotted as a function 
of limiting FIR flux density at 250 /im (solid) , 350 /xm (dashed), 
and 500 fxm (dotted). Detection rates are shown for 3.4 /rm (upper 
left panel), 4.6 fim (upper right), 12 fim (lower left), and 22 fim 
(lower right) sources, where the corresponding 5<r detection limits 
are 0.06, 0.9, and 3.6 mjy, respectively. Note that the large sky 
density of 3.4 fim sources can lead to ambiguity in the identification 
of some sources, so identification rates are typically ~ 10-20% lower 
than the corresponding detection rates. 


(78%) have measured spectroscopic or photometric red- 
shifts < 0.5, suggesting that the low-redshift mode of 
star-forming galaxies dominates sub-mm sources above 
this flux density. Below this FIR flux density, the NIR 
and MIR detection rates drop rapidly as the high-redshiffc 
mode accounts for an increasing fraction of the sub-mm 
sources (see Fig. 2). 

Approximately 1% of the objects in our 250 /nn- 
selected sample (Negrello et al. 2007), including ~ 
50% of sources brighter than 100 mJy at 500 fim 
(Negrello et al. 2010; Hopwood et al. 2011), are expected 
to be strongly-lensed z > 1 galaxies. Of the 46 such 
sources in the G15 region, ~20% are undetected at 12 
and 22 /im, most likely because the 500 /ttn bandpass 
shifts higher in the blackbody curve at high redshift, 
while the rest-frame MIR flux density declines blueward 
of 22 fim (e.g., Rieke et al. 2009). 

The identification of individual NIR and MIR coun- 
terparts is more subtle, as we want to eliminate as 
many of the false matches as possible. Here, we use 
the likelihood ratio technique of Sutherland & Saunders 
(1992), which was implemented for the H- ATLAS sur- 
vey in Smith et al. (2011, hereafter Sll). The likeli- 
hood that a given WISE source is a counterpart to an 
H-ATLAS source is a function of the radial probability 
distribution, /(r), the 3.4 /mi magnitude distribution of 
non-counterparts, n(m), and the 3.4 /im magnitude dis- 
tribution of detected H-ATLAS counterparts, q(m)\ 


q(m)f(r) 

n(m) 

where f(r) is given by 



( 2 ) 


( 3 ) 


Fig. 2. — Distribution of the 3.4-250 /im flux density ratio as a 
function of 250 fim flux density for H-ATLAS sources with WISE 
3.4 /im IDs. Colored points have spectroscopic redshifts from the 
GAMA survey and sources without measured redshifts are indi- 
cated by black points. Dashed lines indicate the approximate lim- 
its of the H-ATLAS and WISE surveys. The majority of sources 
with spectroscopic redshifts (r < 19.4) are well detected by WISE, 
while the no-redshift sources have small flux density ratios that are 
inconsistent with ‘ normal” low-redshift galaxies. 

We measure n(m) directly from the WISE prelim- 
inary release catalog, while we determine q{m) by 
first measuring the magnitude distribution of 3.4 /im 
sources within 10" of the H-ATLAS source positions 
and then subtracting the magnitude distribution of 
background /foreground sources normalized by the area 
within this radius. 

Following Sll, we assume the astrometric uncertainty 
to be circularly symmetric and to have a dependence on 
the 250 fim signal-to-noise ratio (SNR) given by, 


0pos 


0.655 


FWHM 

SNR 


( 4 ) 


where FWHM= 18.1" is. the Herschel beam width. While 
the WISE random astrometric uncertainties are negli- 
gible (~0'/3), the WISE positions in the preliminary 
release catalog are known to suffer from highly non- 
Gaussian systematic offsets that can be as large as 1", so 
we apply a lower limit of cr pos > 1.5". 

Once we have computed the likelihood, Lj , for a pos- 
sible counterpart j, we determine the reliability by sum- 
ming over all possible counterparts, 


Rj — 


Lj 

Li + (1 — Qo) 


(5) 


We consider any 3.4 //m source with R > 0.8 to be a 
reliable counterpart. After performing this procedure 
separately for WISE sources with star- and galaxy-like 
[3.4] — [4.6] colors, we obtain a combined total of 13,898 
WISE counterparts, or 50.6 ± 0.4% of the H-ATLAS cat- 
alog. The expected number of false matches can be ob- 
tained with 

N false ID = £(1 - Ri). (6) 

i 


We estimate that there are 369 false identifications in 
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Fig. 3. — FIR flux density ratio for H- ATLAS sources with WISE 
identifications. Colored grid cells indicate the median redshift of 
a set of modified blackbody models with temperatures uniformly 
distributed between 10 and 60 K and emissivity parameters be- 
tween 0 and 2. Black points indicate sources with 5 <j detections 
in all three SPIRE bands, where open squares have spectroscopic 
redshifts, z < 0.2, and solid triangles have no spectroscopic or 
photometric redshift measurement. Dashed lines indicate the me- 
dian flux density ratios of a stack of sources lacking redshifts. The 
majority of the no-redshift sources appear to be at z ~ 2. Inset: 
Combined distribution of spectroscopic and photometric redshifts 
for H- ATLAS sources with detections at 3.4 (solid histogram), 
12 fim (dashed), and 22 fim (hatched). There are 73 sources at 
z > 1 not shown. 

our sample, corresponding to a contamination rate of 
2.7%. The majority of our IDs have galaxy-like colors, 
but we were able to identify 14 objects with star-like 
colors ([3.4] — [4.6] < 0) that weren’t already identified 
as galaxies with SDSS matching. This corresponds to an 
upper limit of 0.05% for the fraction of H- ATLAS sources 
identified as stars and detected in WISE. 

Comparing the individual WISE identifications to 
SDSS cross-identifications, we find that 84.7% of the 
10, 709 G15 H- ATLAS sources with R > 0.8 SDSS identi- 
fications also had R > 0.8 WISE identifications. The me- 
dian separation between the WISE and SDSS sources is 
076, consistent with the astrometric uncertainties in the 
current version of the WISE catalog. In addition, ~6% 
of these sourcess have separations > 3", consistent with 
the expected contamination rates of the WISE (~ 3%) 
and SDSS (~ 5%) matched samples. 

3. DEMOGRAPHICS 

The redshift distribution of 250 /jm-selected sources 
is likely bimodal, with populations of both moderately- 
star-forming galaxies at 2 < 1 and high- redshift star burst 
galaxies at z ~ 2 (Smith et al. 2011; Lapi et al. 2011). 
Because the latter are typically faint in the observed- 
frame optical and UV, the vast majority of objects with 
measured spectroscopic or photometric redshifts are in 
the low-redshift mode. Of the H- ATLAS sources with 
spectroscopic and photometric redshifts (see the inset of 
Fig. 3), we successfully identify 85.6% in WISE. 

Although most of the WISE identifications with known 
redshifts are at z < 0.8 (>98%), such a sample is 
biased toward low redshift because redshift measure- 


ments require detectable flux in the UV, optical, or near- 
infrared. Approximately one-third (30.4%) of the iden- 
tified sources have reliable spectroscopic redshift and 
two-thirds (66.5%) have photometric redshifts from op- 
tical/NIR data. For the remaining sources, we can esti- 
mate redshifts using their FIR flux density ratios (e.g., 
Amblard et al. 2010; Schulz et al. 2010), which give an 
indication of the location of the peak of the FIR dust 
emission. In Fig. 3, we compare the SPIRE colors of our 
matched catalog to a suite of 10 6 modified blackbody 
spectral energy distributions (Amblard et al 2010). We 
indicate on the plot the median flux density ratios of the 
4, 677 H- ATLAS sources with WISE identifications but 
no redshifts ( ‘no-redshift :5 sample, dashed lines), as well 
as the individual flux density ratios of 311 such sources 
with > 3 a detections in all three SPIRE bands. 

The no-redshift sample has median FIR colors sug- 
gesting z ~ 2, while the subset with three-band SPIRE 
detections may have even higher redshifts, perhaps as 
large as z ^ 3.5. The redshifts derived using the 
Amblard et al. (2010) technique are highly uncertain, as 
redshift is degenerate with dust temperature and it is 
possible that these galaxies have typical dust temper- 
atures different than 35 K (which is typical of high-z 
SMGSj e.g. Chapman et al. 2005). However, in order for 
the no-redshift sample to be at z ~ 0, they would need 
to have very cold dust temperatures (T ~ 12 K) in ad- 
dition to being very faint in the optical (r > 20.5). Al- 
though temperatures as low as 10 K have been seen in 
H- ATLAS galaxies at z < 1 , T ~ 25 - 30 K is more typ- 
ical (Dye et al. 2010). If the no-redshift sample has me- 
dian temperatures typical of these low-redshift H- ATLAS 
galaxies, then our redshift estimate drops to z ~ 1.5. 

Further evidence that the no-redshift sources are at 
z >> 0 can be found in the 3.4-250 /mi ratios (see Fig. 2). 
H-ATLAS sources with GAMA redshifts have NIR-FIR 
flux density ratios that decrease from ~ 0.009 at z ~ 0 to 
~ 0.006 at z ~ 0.35. This trend, which continues toward 
higher redshift, is due primarily to the larger ^-correction 
at 3.5 fim than at 250 /mi (e.g., Rieke et al. 2009). The 
median flux density ratio for the no-redshift sources, by 
comparison, is 0.0023. Although low-redshift galaxies 
will occasionally exhibit such small NIR-FIR flux density 
ratios - Arp 220, for example, has / 3 . 4//250 = 0.0019 - 
such objects are rare, actively star-forming, and unlikely 
to have dust temperatures as small as ^ 15 K. 

We can shed even further light on the properties of 
the no-redshift sample by examining their position in 
WISE color space (see Fig. 4). The technique is sim- 
ilar to those developed for a set of IRAC filters at 
comparable wavelengths (Sajina et al. 2005; Stern et al. 
2005). It was demonstrated that the diagram can be 
divided into stars and early-type galaxies (lower left cor- 
ner, [4.6] - [12] < 1), star- forming galaxies (lower right 
corner), and active galactic nuclei (AGNs), where AGNs 
are selected largely based on their power-law emission 
in the NIR and MIR. This was adapted for WISE by 
Jarrett et al. (2011), 

[4.6] - [12] > 2.2 

[4.6] - [12] < 4.2 
[3.4] - [4.6] < 1.7 (7) 

[3.4] - [4.6] > 0.1 ([4.6] - [12]) + 0.38. 
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Fig. 4. — IR color distribution of ATLAS- WISE matches with 
three-bend WISE photometry, including sources with z < 0.5 
(color coded by redshift) and no-redshift sources (black points). 
The dugrarn can be divided into stars (lower left corner), star- 
forming galaxies (lower right corner), and AGNs (dashed region). 
The color distribution of all three-band WISE detections (including 
sources unmatched to H-ATLAS) is indicated by contours. Sources 
without measured redshifts tend to be redder in [3.4] — [4.6] and 
have a larger AGN fraction than sources with spectroscopic or pho- 
tometric redshift measurements, z < 0.5. 

Of the 4,959 H-ATLAS sources with WISE identifies^ 
tions, r < 0.5, and detections in 3.4, 4.6, and 12 /im, the 
majority are star- forming galaxies, with <0.2% main 
sequence stars (v r hich is consistent with the SDSS/H- 
ATLAS matching done by Thompson et ai. 2010), and a 
small AGN fraction (0.057), 

For the 9% of no-redshift sources with three-band 
WISE photometry, we find a much larger AGN frac- 
tion (0.30), but are subject to a selection bias toward 
AGNs due to their additional flux from warm dust at 
12 /im (Hainline et al. 2011). Furthermore, emission 
from polycyclic aromatic hydrocarbons in high-redshift 
star-forming galaxies can occasionally mimic AGNs in 
their MAC colors, leading to contamination in our AGN 
samples (e.g. Lacy et al. 2004; Donley et al. 2012). For 
comparison, Coppin et al. (2010) performed SED fits on 
a set of blank-field SMG samples and estimated that 
~ 15% of SMGs are dominated by an AGN in the mid- 
infrared. 

4. FIR-NIR PROPERTIES 

Emission from dust in the FIR is often used as an in- 
dicator of the total star formation rate in galaxies (e.g., 
Kennicutt 1998), under the assumption that most of the 
dust heating is provided by young stars and that the 
star-forming regions are optically thick. The 250 /im 
bandpass samples the cooler dust emission, which can 
arise from both star-forming clouds and diffuse regions in 
the ISM (e.g., Eales et al. 2010b; Dunne et al. 2011). By 
contrast, 3.4 /im emission is dominated by stellar sources 
and is a probe of the stellar mass for galaxies without a 
strong AGN contribution (e.g., Stern et al. 2005). 

In Fig. 5, we plot 3.4 /im luminosity as a function of 
250 /ini luminosity for H-ATLAS sources with spectro- 
scopic redshifts, 0.05 < z < 0.3. We k - correct 3.4 /im 



Fig. 5. — Luminosity at 3.4 fim as a function of luminosity at 
250 /im for H-ATLAS sources with spectroscopic redshifts, 0.05 < 
z < 0.3. Both luminosities are rest-frame ^-corrected. Small 
squares are colored according to their redshift and red stars in- 
dicate AGN, selected with the criteria given in Eq. 7. Dashed lines 
indicate the approximate selection limits for sources at z — 0.05 
and 0.2. The majority of the AGN lie above the relationship, sug- 
gesting that L 3.4 contains a non-negligible contribution from hot 
dust emitting near the central black hole. 

using a power law interpolation of the [3.4] — [4.6] color. 
For the 250 /im flux densities, the majority of the galax- 
ies in our sample do not have sufficient FIR flux den- 
sity to obtain reliable dust temperatures, as is needed 
for a proper k correction. Instead, we use the median 
dust temperature of 26 K found by Dye et al. (2010) us- 
ing PACS and SPIRE data from the H-ATLAS SDP of 
1,346 0.1 < 2 < 0.5 sources. Although the r < 19.4 limit 
of the GAMA spectroscopic survey restricts the range 
of luminosities that can appear on the plot, we estimate 
that GAMA misses only ~ 5% of matched z < 0.3 galax- 
ies based upon the deeper (r < 20.5) subsample of H- 
ATLAS galaxies v.ith photometric redshifts. 

We find that the majority of the IR-color-selected AGN 
(red points) lie above the relationship for star-forming 
galaxies, as L3 4 will contain a non-negligible contribu- 
tion from hot dust emitting near the central black hole 
in addition to emission from low-mass stars. However, if 
we exclude AGNs, we find an approximately linear cor- 
relation, with a best-fit power law index, a = 0.98io]os- 
The quoted systematic uncertainties on the power law 
index were determined by allowing for a range of possi- 
ble dust temperatures, 18 < T < 34 K, when performing 
250 /im ^-corrections. The intrinsic scatter about this 
relationship is 0.18 ± 0.01 dex, or ~ 50%). 

The existence of a linear correlation between L250 and 
L3.4 across one and a half decades in luminosity suggests 
a close relationship between the cold dust probed by the 
FIR and the stellar mass probed by the NIR. Previous 
indications with Herschel have shown the 250 fi m lumi- 
nosity density to tightly correlate with both the 24 /im 
luminosity density, a star formation rate indicator, and 
the total infrared luminosity (Elbaz et al. 2010). This 
fact, coupled with a linear correlation between stellar 
mass and star formation rate for the general star-forming 
galaxy population (Daddi et al. 2007; Elbaz et al. 2007; 
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Noeske et al. 2007, Donoso et al., in prep), suggests that 
the majority of the low-redshift H-ATLAS sources are 
actually “normal” star-forming galaxies. 
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and the Jet Propulsion Laboratory/ California Institute 
of Technology, funded by the National Aeronautics and 


Space Administration. 

Herschel is an ESA space observatory with science in- 
struments provided by European-led Principal Investiga- 
tor consortia with sig- nificant participation from NASA. 
U.S. participants in Herschel ATLAS acknowledge sup- 
port provided by NASA through a contract issued from 
JPL. 

In addition, we thank Dan Stern and Roberto Assef 
for helpful discussions about the quasar selection. 


REFERENCES 


Abazajian, K. N. et al. 2009, ApJS, 182, 543 
Amblard, A. et al. 2010, A&A, 518, L9+ 

Chapman, S. C., Blain, A. W., Smail, L, & Ivison, R. J. 2005, 
ApJ, 622, 772 

Coppin, K. et al. 2010, ApJ, 713, 503 

Daddi, E. et al. 2007, ApJ, 670, 156 

Donley, J. L. et al. 2012, arXiv:astro-ph/ 120 1.3899 

Driver, S. P. et al. 2011, MNRAS, 413, 971 

Dunne, L. et al, 2011, MNRAS, 417, 1510 

Dye, S. et al. 2010, A&A, 518, L10+ 

Eales, S. et al. 2009, ApJ, 707, 1779 

— . 2010a, PASP, 122, 499 

Eales, S. A. et al. 2010b, A&A, 518, L62+ 

Elbaz, D. et al, 2007, A&A, 468, 33 
— 2010, A&A, 518, L29+ 

Griffin, M. J. et al. 2010, A&A, 518, L3+ 

Hainline, L. J., Blain, A, W., Smail, L, Alexander, D. M., Armus, 
L., Chapman, S. C., & Ivison, R. J. 2011, ApJ, 740, 96 
Hopwood, R. et al. 2011, ApJ, 728, L4+ 

Jarrett, T. H. et al. 2011, ApJ, 735, 112 
Kennicutt, Jr., R. C. 1998, ApJ, 498, 541 


Lacy, M. et al. 2004, ApJS, 154, 166 
Lapi, A. et al. 2011, ApJ, 742, 24 
Lawrence, A. et al. 2007, MNRAS, 379, 1599 
Negrello, M. et al. 2007, MNRAS, 377, 1557 
— 2010, Science, 330, 800 
Noeske, K. G. et al. 2007, ApJ, 660, L43 
Pascale, E. et al. 2008, ApJ, 681, 400 
Pilbratt, G. L. et al. 2010, A&A, 518, L1+ 

Poglitsch, A. et al. 2010, A&A, 518, L2+ 

Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., Perez-Gonzalez, 
P. G., Blaylock, M., Donley, J. L., & Maxcillac, D. 2009, ApJ, 
692, 556 

Rigby, E. E. et al. 2011, MNRAS, 415, 2336 
Sajina, A., Lacy, M., & Scott, D. 2005, ApJ, 621, 256 
Schulz, B. et al. 2010, A&A, 518, L32+ 

Smith, D. J. B. et al. 2011, MNRAS, 416, 857 

Spergel, D. N. et al. 2007, ApJS, 170, 377 

Stern, D. et al. 2005, ApJ, 631, 163 

Sutherland, W. & Saunders, W. 1992, MNRAS, 259, 413 

Thompson, M. A. et al. 2010, A&A, 518, L134+ 

Wright, E. L. et al. 2010, AJ, 140, 1868 



